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Abstract 

We show that the probabihty mass function of the riff-shuffle distri- 
bution, also known as the minimum negative binomial distribution, is 
unimodal, but in general not log-concave. 

Keywords: Discrete distribution, unimodal 

2000 Mathematics Subject Classifications: Primary: 60C05; Secondary: 11K31. 



The rifF-shuffle distribution is defined for < p < 1 and integers m > 1 by 
the probability mass function 

fk = ' ^) + g V), fc = 0, . . . , m - 1, 

where p + q = 1. If cards are taken with probabilities p and q from two decks 
of m cards, then the number of cards that have been taken from the remaining 
deck when one deck has been depleted follows this distribution W . Equivalently, 
it is the distribution of the random variable X — m^ where X denotes the number 
of Bernoulli trials with individual success probability p until either m successes 
or m failures have been observed. This explains why it is also called minimum 
negative binomial distribution [3.J. Johnson, Kotz, and Kemp [1] report that 
Uppuluri and Blot's [2] plots exhibit unimodality of the distribution. In this 
short note we confirm this observation. 

Recall that a sequence (5^)^=0 called log-concave if gf. > gu-igk+i holds 
for 1 < fc < n. It is well known that a positive log-concave sequence is unimodal. 
The riff-shuffle distribution, however, provides examples of unimodal sequences 
that are not log-concave. 
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Proposition 1. The probability mass function of the riff-shuffle distribution is 
log-concave for p = \ . For p ^ ^ and m sufficiently large, it is not log-concave. 

Proof. After cancelling the binomial coefBcients, using the relation 

u\ u fu — 1 
V J V \v — 1 

the log-concavity condition > fk-ifk+i becomes 
For p= \, this reduces to 

k{m + k) 



1 > 



(A;+l)(m + A;-l)' 



which is easily seen to be true for fc = l,...,m— 1. To show the second assertion, 
we set A: = m — 1. After some simplifications, we see that ^ /m— 2/m Is 

equivalent to 

2(2m - l)(m - l){pqf'^-^ + (m - l)(p2™g2m-2 ^ 2{pqf"'-^ + p2"-2^2m) 

>(2m-l)(m-l) (^^ + ^)(pg)2— 1. 

The left-hand side of this inequality is ~ Ann? {pq)'^"^~^ as m — » 00, whereas 
the right-hand side is ~ 2{p/q + q/p)m'^{pq)'^"^~^ . Since p/q + q/p is a convex 
function of < p < 1 with minimum 2 at p = i, the right-hand side is larger 
than the left-hand side for p ^ ^ and sufficiently largo m. □ 

Since the sufficient criterion of log-concavity fails us for p 7^ ^ , we now prove 
directly that the distribution is unimodal. 

Proposition 2. The probability mass function of the riff-shuffle distribution is 
unimodal, with the maximum max fk occuring for at most two (adjacent) values 
ofk. 

Proof. For p = 5 it is straightforward to show that /o < /i < • • • < fm-2 = 
fm-i- Now we assume w.l.o.g. that Q < p < i. Defining the auxiliary function 
h{x) := p™q^ + q™p^, we have 

fk+i _ k + m h{k + l) 
fk - k + 1 h{k) ' k-0,...,m-2. 

Therefore, fk+i < fk is equivalent to g{k) > 0, where 

h{x) X + m 



g{x) :-- 



h{x + l) x + l 
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We show that g is concave on [0,m — 2]. The second derivative 

d^ h{x) _{q- p)(logp - logg)2p™+^g™+^(p'"g^+i - p^+ig'") 

h{x + l) ~ (^pmqX+l j^px+lqm^3 

is non-positive for < x < m — 2, since {p/q)™^^^^ < 1 implies p™q^+'^ — 
px+iqun ^ g_ xhus, h{x)/h{x+ 1) is concave. Since {x + m)/{x+ 1) is a convex 
function of a; G [0, m — 2], the function g is indeed concave. At a; = m — 2, the 
function value 

h(m-2) 2m -2 m-2 ^ m m-2 

o(rn - 2) = — — = = ^ - 2 

^ ^ h{m-l) m-l pmqm-l ^ qmpm-l 

is positive. From this and concavity we conclude that, if g{l) > for some I, 
then the values g{k), I < k < m — 2, are positive. In other words, as soon as 
there is a descent < the sequence (//c)fe!^+i strictly decreases. □ 
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